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We present an alternative scenario for cosmic structure formation where initial fluctua- 
tions are due to Kalb-Ramond axions produced during a pre-big bang phase of inflation. We 
investigate whether this scenario, where the fluctuations are induced by seeds and therefore 
are of isocurvature nature, can be brought in agreement with present observations by a 
suitable choice of cosmological parameters. We also discuss several observational signatures 
. which can distinguish axion seeds from standard inflationary models. We finally discuss 

' the gravitational wave background induced in this model and we show that it may be well 

within the range of future observations. 

(D ■ PACS Numbers : 98.80. Cq 98.80.E 

^ ; I. INTRODUCTION 

> 

It is commonly assumed that an inflationary phase is necessary in order to construct a consistent cosmo- 
logical model. The familiar adiabatic inflationary scenario owes its popularity to the fact that it solves the 
horizon and flatness problem and at the same time provides a consistent model for the origin of cosmological 
' perturbations. In particular, it naturally leads to a flat (Harrison-Zel'dovich) spectrum of perturbations on 

I large scale and to coherent acoustic oscillations on intermediate scales which manifest themselves as "peaks" 

■ in the Cosmic Microwave Background (CMB) anisotropies. 

After the recent measurements of the intermediate scale CMB anisotropy power spectrum (Hj^-^, flat 
Q^. adiabatic models seem to be favored Nevertheless none of the many inflationary scenarios which 

have been developed during the last 20 years has been constructed consistently on the bases of a serious 
theory of high energy physics; inflation has always been seen as an effective model pointing to a greater 
more fundamental theory which has not been clarified so far. We believe that superstrings are presently the 
^ , most promising candidate for such a theory but on the other hand it is well known that it is not possible to 

derive an inflationary model from a string theory effective action on a generic background, the reason being 
that the non-minimal coupling between the dilaton and the metric slows down the expansion of the universe 
' spoiling the solution of the problems for which inflation has been invoked. 

H I The pre-big bang idea [§| represents in this context one of the first and most interesting attempts to develop 

a new cosmological scenario which solves the horizon and flatness problems, based on string theory. In this 
radically new picture, the underlying duality symmetry pO| | present in the low energy sector of string theory 
naturally selects perturbative initial conditions and automatically leads to an inflationary phase prior to the 
big bang during which curvature and the dilaton are growing Besides its many appealing features, 

this scenario is known to face several problems such as the lack of a complete and consistent description of 
the high coupling and high curvature regime where the transition between the pre-big bang and the post-big 
bang phase and the stabilization of the dilaton should take place fl^ j. Furthermore, opinions vary as to 
whether the initial conditions in the pre-big bang need a large amount of fine tuning ||Tl] , [l3|| . On a more 
phenomenological side, it is nevertheless important to study whether this scenario can provide the features 
that we observe in the universe today. For recent review articles discussing several of the previous points we 
refer the reader to |14|Jl5[| . For a comparison of the pre-big bang model with new cosmological models based 
on string theory see ||lq| T 

A realistic cosmological model has to generate large-scale matter perturbations and to reproduce the slope 
and the amplitude of CMB anisotropy spectrum. The pre-big bang scenario was thought for some time to be 
unable to provide a scale-invariant spectrum of perturbations. First-order tensor and scalar perturbations 
in the metric, as well as perturbations of the moduli fields, were found to be characterized by extremely 
"blue" spectra jl^. This large tilt, together with a natural normalization imposed by the string cutoff at 
the shortest amplified scales, make their contribution to large-scale structure completely negligible. 
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However, it was later realized that the spectral tilt of the axion, a universal field in string theory, can 
assume a whole range of values depending on the behavior of the internal and external dimensions and in 
particular it can naturally provide a scale- invariant spectrum of perturbations jl^ This result reopened 
the possibility that pre-big bang cosmology may contain a natural mechanism for generating large-scale 
CMB anisotropics via the "seed" mechanism 

This possibility was analyzed in Refs. ||2^ , |23|| for massless axions and in Ref. |2j] for very light axions. 
These analytical treatments are restricted to large angular scales. We have extended the study to smaller 
scales with the help of numerical calculations. First results for this work have been reported in a letter 
p5[ , where a strong correlation between the axion spectrum, 71^, and the height of the peak was noticed. 
A range of values around n„ = 1.4 (slightly blue spectra) appeared to be favored by a simultaneous fit to 
the normalization on large angular scales observed by COBE |26| and the data on the first acoustic peak 
available at that time. 

In this companion paper we present a full explanation of the details of these calculations for the CMB 
angular power spectrum and for the dark matter power spectrum and we study the problem of the "de- 
coherence" of axion perturbations which has been ignored in the previous work. Furthermore, we expand 
on the results on the observational signatures presented in the Letter |2^] and we discuss them in the light 
of the new CMB anisotropy data presently available by investigating the cosmological parameter-space of 
the model. We also discuss CMB polarization for our model and the contribution of the gravitational wave 
background induced by axion perturbations. 

We will suppose for the rest of this work that modes with frequencies relevant to CMB physics are 
unaffected by the the transition from the pre- to the post-big bang phase. This may seem a strong assumption 
given our ignorance on the graceful exit and on the duration of the intermediate string phase. Indeed, 
the nature and features of the transition are still unknown and may lead to important changes on the 
characteristics of the axions spectrum, but we are confident that this can only affect scales much smaller 
than those we are interested in: For this to hold, we just need that, at the moment we enter the string phase, 
the CMB modes are well outside the horizon and therefore not altered by any causal process. As shown 
in p?! using general arguments, the long-wavelength part of the solution to the equation of motion of the 
axion field is always dominated by the contribution of the frozen modes even if the background evolution 
includes a high-curvature phase in which the action and the perturbation equations are not known. This 
argument has been tested by numerical examples in |28| for a gravitational wave background spectrum. 
Recently a graceful exit model considering general high-curvature and coupling corrections has appeared in 
the literature p9t and numerical tests on the axion spectrum have been performed and will be published 
soon pO{ to confirm our assumption. 

The paper is organized as follows: in the next section we study axion production in the pre-big bang 
and explain the details of the computation of the axion energy-momentum tensor which plays the role of 
the "seed" in our model. In Section III we determine the CMB anisotropy and dark matter spectra. We 
study the problem of decoherence and show that the coherent approximation is very good for this model. 
In Section IV we compare our result with CMB and Supernova data and present a cosmological parameter 
estimation for this scenario. We also examine and discuss the normalization and the kink in the axion 
spectrum which is required to fit observations. Section V is devoted to a novel prediction of axion seeds: 
the tensor component of their energy-momentum tensor induces a gravity wave background which might be 
observable. In Section VI wc summarize our conclusions. 



II. AXION SEEDS FROM STRING COSMOLOGY 



A. Extra dimensions in string cosmology 



The minimal low energy effective action of the NS-NS sector of string theory in the string frame is given 



by [0 



S 
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i?io + (V0io)' - -Hi 



(2.1) 



where we have included the 10-dimensional antisymmetric tensor H^^a = d[i_iB^a], but no gauge or fermion 
fields. 

We assume that the 10-dimensional metric can be factorized into a "large" four-dimensional part and a 
"small" six-dimensional metric, 

ds^ = g^^dx^'dx" + e^'^SijdX'dX-' (2.2) 

(/i, = 0, . . . , 3 and I, J = 1, . . . , 6), where /3 depends only on time, (3 = P{t). If the six dimensional piece is 
compactified to a very small radius, the lowest energy Kaluza-Klein modes yield the four-dimensional action 



S ^ d xy'\g\e 



R + (V0)2 - 3(V/3)2 - ie2'*(VCT)^ 



(2.3) 



Here we have introduced the four-dimensional axion field a defined by 

jjt^va ^ e'^e'"'"''V;3cr. (2.4) 



The action (2.3) and the definition (2.4) include the four-dimensional dilaton field, (/>, the pseudo-scalar axion 
field, a, which represents the degrees of freedom of the antisymmetric three tensor field H, and a modulus 
field, /3, which parameterizes the radius, or the "breathing mode" , of the six-dimcnsional internal space. 
Like the dilaton, also the axion field (not to be confused with the Peccci-Quinn axion) is universal in string 
theory. 

Let us assume a homogeneous dilaton background, (j) = (j){t)^ and an external four-dimensional spacetime 
adequately described by a standard, spatially flat FLRW metric with scale factor a{t), 

g^, = diag[-l, a\t), a^{t),a\t)]. (2.5) 

In the following we shall also make use of the metric 

g^, =a2(77)diag[-l,l,l,l], (2.6) 

where we have introduced the conformal time "q given hy drj — dt/ a (we shall use an over-dot to indicate 
a derivative with respect to conformal time, ' = d/drj). With our choice of the external metric, the four- 
dimensional dilaton is related to the 10-dimensional one by 

= 010 - 6/3. (2.7) 

When the axion field is trivial, it = 0, or its contribution to the global dynamics of the universe is negligible, 
the equations derived from the action (pTsh are invariant under duality transformations. 



a{t) l/a{-t), (j){t) (j){^t) - 61n(a(-i)). (2.8) 

This invariance {scale factor duality) represents one of the key motivations behind the pre-big bang sce- 
nario [^. The field equations for a,(j) and (3 are solved by the following power laws, known as dilaton- 
vacuum solutions in the pre-big bang for rj < —rji: 



"C'n^J . , .-^^^J , (2,9) 

where 5 and C satisfy the Kasner constraint, 

3(5^+6C^ = l. (2.10) 

Here —rji is the (conformal) time at which curvature and dilaton become so large that loop corrections from 
string theory have to be taken into account. It is hoped that these corrections then lead to a radiation 
dominated Friedman universe with "frozen" dilaton aX r] > rji. From these solutions one can see that, during 
the pre-big bang phase, i.e. for negative conformal time 77, a negative 5 and a positive C are required to make 
the external three-dimensional space expand and the internal six-dimensional space contract. Therefore 5 
has to lie in the interval — 1/\/3 < (5 < 0, which leads always to a growing dilaton and growing four-curvature, 
R - (a/a2)2 oc l/{a-qf cx (-77)^^. 
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B. Amplification of axion quantum fluctuations 



In this subsection we briefly review the mechanism for the generation of a primordial quasi-scale-invariant 
spectrum from the pre-big bang phase and we discuss the dependence of the spectral index on the evolution 
of the internal and external dimensions of the pre-big bang universe. Using as initial conditions the axion 
field obtained during the pre-big bang phase, we then analyze its evolution after the big bang in a critical 
FLRW universe with and without cosmological constant, paying particular attention to the frequency modes 
that enter into the calculation of the CMB anisotropy power spectrum. 

As in previous works l22| - [25| we suppose that the contribution of the axion field to the equations of 
motion for (j), a and (3 is negligible and that the evolution of the dilaton, the moduli, and the scale factor are 



governed by the dilaton- vacuum solutions (2.9). Nevertheless, quantum fluctuations of all the fields are of 
course present and we will show that quantum fluctuations of the axion field can seed density perturbations 
and CMB anisotropics in the post-big bang era. To this goal we have to study the axion evolution equation 
and the spectrum of axions produced during the pre-big bang phase due to their coupling to the background 
gravitational field and the dilaton. 



Varying the action (|2.3| ) with respect to the field a in the string frame yields the equation of motion 

V^(e"^V''cr) = 0. (2.11) 

The study of this equation is conveniently performed by using the canonical variable given by 

= a^cr = ae'^/V, (2.12) 

which "diagonalizes" the perturbed action expanded up to second order. The factor aA is the so called pump 
field of the axion. The Fourier modes 'ijjk{ri) satisfy a canonical linear second-order equation, completely 
decoupled from the other fields, 

4 + (^fc' - ^) V^k = 0. (2.13) 
This is the evolution equation for the axion field. 



Eq. (2.13) is equivalent to the equation for a classical harmonic oscillator with parametric evolution driven 
by the time dependent mass term oaI aA- When the time evolution of the velocity of the pump field, a a, is 
sufficiently slow such that, for a given mode fc, OAjaA <C fc'^, we are in the adiabatic regime with the result 
that no particles are created. When the acceleration in the pump field is high enough to violate the adiabatic 
regime, quantum particle production starts. The evolution of the axion field and the resulting spectrum of 
particles are fully determined by the time behavior of the pump field in the different phases of the universe. 
In particular, a strong difference in this behavior exists between the pre-big bang phase and the standard 
radiation and matter dominated eras in the post-big bang universe. 

The pre-big bang phase is characterized by an accelerated evolution of the pump field, 

a^(x(-,y)^ (2.14) 



wher e 8 < is the power which characterizes the evolution of the external dimensions, Eq. (2.9). Using 



Eq. (2.13), the evolution equation of the axion can be written as 



V-k + fcM 1 - ^^^1^ ) = 0, (2.15) 



where x = krj. This equation is solved analytically in terms of the Hankel functions rf^^'^H^^^ and rj^/'^Il'if' 
with ^ = |7 - 1/2|. 

At very early time, a per turbation of given wave number k is well inside the horizon, \x\ = \kri\ 3> 1, and 



the solutions of Eq. (2.15) are harmonic oscillations which can be consistently normalized to the vacuum 



fluctuation spectrum for 77 — s- —00. This initial condition implies that the H^i^ mode is absent and 
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Mv) = {-vY'^H^\k'n), M = ^ - 7 - for 77 < -r;i. (2.16) 

At 77 = — 771 wc "glue" the pre-big bang phase to a radiation dominated post-big bang era starting at 77 = 771. 

After the singularity, during the standard radiation and matter dominated eras, the dilaton is frozen, 
(j) = const, and the pump field is proportional to the standard scale factor, aA oc a. The scale factor, a, 
and its second derivative, a, are given by Friedman's equations. For a critical universe, which we consider 
throughout our calculations and which is certainly a good approximation until redshifts 2: < 5, we have 

a AttG , 2a^A 

a\p-3p) + ^, 2.17 

SttG o a^A 

- = — a2p+— . 2.18 
3 3 

Energy conservation for radiation (r) and matter („) yields p,- cx 1/a^ and p,„ oc 1/a'^, with p = + Pm 
and p = pr/i] Pr is the radiation energy-density, pm is the matter energy-density, and p the pressure of the 
radiation fluid. At early times, when A is negligible, these equations have a simple analytical solution, 

I I \2\ ^ ( 3 V''^ 77e 



where 772^ is the transition time between the radiation and the matter dominated era, Pr{ileq) = Pmi^eq) = 
Peql"^- The mass term during the post-big bang becomes 

CLA d \ , . 

— = - = o (2.20) 

When A is non vanishing, the solution for the mass term can be found numerically but since the contribution 
of a small cosmological constant to the scale factor becomes important only at late time, the solutions to 



(2.15) are almost unaffected; this has been checked by numerical tests. In the radiation dominated era, 
77 < ?7eq, the mass term can be approximated by a a/ a a — 1/(277*7;). 

We now study the axion evolution in the post-big bang era. Let us write the term in parenthesis on the 



left hand side of the axion equation of motion, Eq. (2.13|), as 



a 



In order to study the solution of Eq. (2.13) we have to study the ratio of the dimensionless mass term 
{d/a)r]'^ and to be compared with unity. As long as we are well in the radiation dominated era, rj <C 7/*, 
the dimensio nless mass term is small and particle creation induced by the pump field is negligible at early 
times. Eq. ( ^.13 ) then is a harmonic equation solved by free plane waves. 



MV) = ^[c+{k)e-'''' + c_(k)e'^'']. (2.22) 



By matching the two solutions (|2.16| ) and (|2.22| ) at the transition time 771 we obtain, for |fc77i| ^ 1 and 
7?eg » 77 > r]i, 



so that 




c±(k)=±c(k), with (|c(k)|^)~<; yklj ^^^^ (2.23) 

k > ki , 



c(k) 

V^k = -^sin(A:(77 - 77i)), for 771 < 77 < ?7eg . (2.24) 

Here fci — I/771 represents the maximal amplified frequency of the pre-big bang phase. As already discussed 
in the introduction we suppose that modes with frequencies much lower than ki are unaffected by the 
unknown details of the transition from the pre- to the post-big bang phase. 
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The energy-density distribution of the produced axions is then 



The axion spectral index n^- is related to the power which characterizes the evolution of the external dimen- 
sions by 

71^ =4-2^ = 3 + 27 = 2 , (2.26) 



which follows from Eq. ( 2.16 ). In order not to over-produce infrared axions we have to require fi < 3/2, 
or ricr > 1, which implies S > —1/3. As already pointed out in p2| , the limiting value /i = 3/2 corresponds 
precisely to a Harrison-Zel'dovich spectrum of CMB anisotropics on large scale. In terms of the evolution of 
the scale factor, this corresponds to an isotropic expansion and contraction respectively of the external and 
internal dimensions, 

a(x ^ (X (-7-1)-'^/^. (2.27) 

Notice that only for a 10-dimensional spacetime, symmetrical expansion and contraction corresponds to a 
flat axion spectrum which induces a Harrison-Zel'dovich spectrum of CMB fluctuations ! 

Nevertheless, as will be discussed in Section IV, at very large scales and very early (negative) times, we 
will need a slightly blue axion spectrum to fit CMB data. This requires a somewhat larger value of S, i.e. a 
slower expansion of the external dimensions and, correspondingly, a somewhat faster contraction of internal 
dimensions at early time. This blue spectrum cannot be maintained up to the string scale because the fixed 
normalization at the string scale to gf = [(fci/ai)/r7ipianck]^ ~ 0.01 + 10~^ would lead to much too small 
amplitudes at the COBE scale. 

Let us therefore investigate what happens if the universe expands with some expansion law described by 
6- at early times, t] < rjf, > —rji and then switches to an expansion law given by i5+ after Sufficiently 
short wavelength modes which are inside the horizon during the entire epoch rj < rjb, which satisfy \krif,\ < 1, 
are not influenced by this change in the expansion law. The term Ha/cla is indeed sub-dominant in the 



equation of motion for t/^k during this epoch and hence the B ogoliubov coefhcicnt |c(k)p of Eq. (2.23) is not 



influenced by the transition; we just obtain the result (2.23) with /i = 



The situation is different if a mode exits the horizon before rji,. Then the "incoming" solution '0('/ < 
77;,) = {—r])^/'^Hj^}{kri) differs from the vacuum solution and matching it to the general "outgoing" solu- 

tion, ^Pir] > r^b) = h{-r,y/^Hj,'^^{kr,) + b2{-v)^/^Hj^}ikr]), yields 62 - 61 = (r(^_)/r(^+)) \kr,b/2\''+-^- . 
Correspondingly, the coefficient |c(k)p is changed by a factor I62 ^ In a model where the expansion 
law changes at a well defined time rji, = — l/fcfc, we therefore get the following Bogoliubov coefficients in the 
post-big bang radiation era (see Fig. IT]) : 



k\ 1 J (^k/k^)2^+-2^- for < fc^ 



c(k)n-^^j ^r-' -l^l (2.28) 

We do not want to specify the event which may have triggered such a transition from ^^(fc < k(,) = 
4 — 2/i_ = 1 + e to n^ik > kf,) — 1, but there are certainly different possibilities. For example, it is 
interesting to note that isotropic expansion and contraction, a oc 1/6, in a 26-dimensional spacetime gives 
S- — 1/5, or Ua = 1.33, which corresponds to £ = 1/3, just about the "tilt" needed to fit the observed 
CMB anisotropics (see Section III). Therefore, if we start out the pre-big bang phase with a 26-dimensional 
bosonic string vacuum (which we know to be unstable due to the presence of tachyons) which then "decays" 
to a supersymmetric and 10-dimensional string vacuum at some time r]i,, which corresponds to a comoving 
energy scale k},, this could induce the required tilt. 

We now study the modification in the axion spectrum during the post-big bang era, where = a. As 
we have seen above, during the radiation era, rj < r]^, the dimensionless mass term is small. Furthermore, 
once a mode enters the horizon, kij > 1, the fc^-term always dominates over the mass term and there is no 
more particle creation. Therefore modes which enter the horizon before equality, krj^ > 1, are not amplified 
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LogfkJ Log[k] 

FIG. 1. Evolution of the axion spectral index during the pre-big bang. The value gi is the string coupling constant 
given by the string scale divided by the Planck scale, g\ — (fci/ai)/mpianck (see Section IV). 



any further in the post-big bang phase. The spectrum of axion perturbations for these modes remains 
unaffected. However, the low frequency tail of the spectrum is further modified as soon as we enter in the 
matter dominated era, where the dimensionless mass becomes of order unity. The modes which enter the 
horizon after equahty, fcry* < 1, are amplified. This amphfication of low frequency modes has important 



consequences on the angular spectrum of the CMB as will shall discuss in detail in Subsection III B 

The behavior of the dimensionless mass term (0/0)77^ together with two modes that enter the horizon 
before and after equality have been plotted in Fig. H. As one can see, only modes entering the horizon 



after equality are amplified. 
Eq. (2.13) becomes 



Deep in the matter era 77 ^ rjeq, the dimensionless mass term is constant and 



fc' - ^ ) = 0. 



This equation can again be solved in terms of Hankcl functions, 

Mv) = v'^^[AH^%ikr,) + BH^]l{k7j)], for 77 » 77^,, 



(2.29) 



(2.30) 



where A and B are cons tants to be determined by matching conditions (see (2^). The post-big bang 
solutions ( 2.24 ) and ( 2.30 ) are only correct far from matter-radiation equality ijeq and in order to compute 
CMB anisotropics we require better precision for these solutions also for 77 ^ rjeq. We therefore solve the 
axion equation of motion Eq. (2.13) numerically, from the early radiation era through the radiation-matter 
transition. 

The axion field is then given by 



cr(k,77) 



V'k(?7) -^^"^(fc,??), 



i\/k 



where the variable ip is the solution of equation 



(p+ fc2--U = 



with initial condition (obtained from the pre-big bang solution Eq. (2.24)) 

(p(fc, 77) = sin(fc77), 77 <C77». 



(2.31) 



(2.32) 



(2.33) 
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FIG. 2. The dimensionless mass term, {a/ajr/^ , (thick hne) and two modes that enter the horizon before and after 
equality. The mode which enters the horizon before equality, tpik = 10fceq,r?), is unaffected by the pump field and 
begins to oscillate without being amplified. The mode that enters the horizon after equality, tp[k = 0.1fceq,r;), is 
amplified and begins to oscillate later. 



We have solved Eq. ( 2.32 ) numerically in this work using the mass term ( 2.2C ). The pre- factor c(k) is a 
stochastic Gaussian field with power spectrum 



(|c(k)p) = (k/hY 



(2.34) 



where n„ is again the primordial spectral index (2.26), our free parameter which depends on the higher- 
dimensional pre-big bang phase. 



C. Axion quantum fluctuations as seeds 

We are now ready to consider the axion field as a source of the linear cosmological perturbation equations. 
As in previous works [p2|-p5[ we suppose that the contribution of the axions to the cosmic fluid can be 
neglected and that they interact with it only gravitationally. They then play the role of seeds which, 
by their gravitational field, induce fiuctuations in the cosmic fiuid pl[ |. The back-reaction of the metric 
perturbations on the evolution of seeds is second order and can be neglected in first order perturbation 
theory. The evolution of axions can be computed by using the solutions of the axion field equation in the 



unperturbed background geometry, Eq. (2.13) 



The axion field cr is a Gaussian stochastic variable. Its contribution to the perturbation equations is given 
in terms of its energy-momentum tensor, 

T^^) = d^ad.a - ^g^AdaCrf, (2.35) 

which is quadratic in a and therefore not Gaussian. Moreover, although the axion field evolves according to 
a linear equation, it will enter into the perturbation equations through T^^-* which evolves non-linearly. 

The perturbations in the dark matter and radiation components are set to zero in the initial conditions 
and are subsequently induced by the gravitational field of the axion. Hence, axion seed perturbations belong 
to the class of isocurvature perturbations. However, they differ from topological defects by being "acausal", 
i.e. they have non- vanishing correlations on super-Hubble scales, since they are due to field excitations 
induced during an inflationary era. 

As we have seen above, the axion power spectrum obeys a simple power law with cutoff and is in general 
not analytic at A; = 0. Furthermore, axion perturbations do not, in general, display the scaling behavior 
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expected from topological defects. In the pre-big bang we have an additional scale, the string scale ki, which 
breaks scale-invariance. The axion spectrum on large scales is therefore not determined by dimensional 
arguments since there are dimensionless factors of the form (fc/fci)" which may alter the spectrumQ The 
significance of these points will become clearer later in the paper. 

As in p5[ |, we first consider a critical universe (total density parameter = 1) consisting of cold 
dark matter, baryons, photons, and three types of massless neutrino, with or without a cosmological 
constant. We choose the baryonic density parameter flf, — 0.05 and the value of the Hubble parameter 
Ho = 100/ikms-iMpc^^ with h = 0.65. 

The linear perturbation equations for this universe in Fourier space are of the form 

VX^S, (2.36) 

where X is a long vector containing all the fiuid perturbation variables which depends on the wave number 
k and conformal time ?], 5 is a source vector which consists of certain combinations of the seed energy 
momentum tensor an d V is a linear ordinary differential operator. More details on the linear system of 
differential equations ( 2.36 ) can be found in js^ and references therein. 

For a given initial condition, this equation can in general be solved by means of a Green's function, g{ri, if), 
in the form 

X{k,7jo)^ e(k, 770, 77)5(k, 77)^77. (2.37) 

We want to determine power spectra or, more generally, quadratic expectation values of the form 

(X,(k,r;o)X,(k,77o)*), (2.38) 



which, according to Eq. (2.37), are given by 

(X,(k,77o)X,(k,77o)*) = / / g.aivo,ri)g*^{vo,ri'){Siirj)S*„Md7^d7j'. (2.39) 



Vir, "'') 



(Sums over double indices are understood.) 

We therefore have to compute the unequal time correlators, {Si{ri)S*-^{r]')), of the seed energy-momentum 
tensor. This problem can, in general, be solved by an eigenvector expansion method |32 3^, as it will be 
done in Subsection III E| . However, if the source evolution is linear, the problem becomes especially simple. 
In this "coherent" case, we have 

^jiv) = Fji{ri,riin)Si{r]in), (2.40) 

with a deterministic transfer function Fji. In this situation we can, by a simple change of variables, diago- 
nalize the hermitian, positive initial equal time correlation matrix, 

{Sl{riin)Sm{'ntn)) = h^lm- 

Inserting this in Eq. ( 2.39| ) yields 

{X,{r])X;i7j')) = (^J\u{m,v)F^l{r],v^n)^/Xldr^ G jm{ilQ,v')F,M ,r^,n)y/K.dri'^ Si^. (2.41) 



We therefore obtain exactly the same result as the one obtained by replacing the stochastic variable S by 



the deterministic source term 5''''^*) given by 



Actually the radiation - matter transition scale 7;* represents a scale which is also present in models with topological 
defects, but deep in the radiation or matter era this scale has no significance, whereas as we shall see the above factors 
multiply the entire power spectrum of fiuctuations. 
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sl'''\v)S\''''\rj) = F,,{rj,n,^)Fa{v,mn)\i = exp{te,,)^ (\SM\')mvW), (2.42) 

where 9ji is a, in principle unknown, phase which has to be determined case by case. Clearly 9jj — 0. When 
the stochastic variable iS is real (as in our case) exp(idji) = ±1. This linear or coherent approximation will 



be fully used in this paper. We shall test its validity in Subsection [II B . 

It is useful to split the energy-momentum tensor of the axion seeds (2.35) into a scalar, vector, and tensor 
part since the perturbations generated by each of these components evolves independently. Due to statistical 
isotropy these three modes are uncorrelated. This also corresponds to a decomposition of the source term 
S into a scalar, vector, and tensor contributions, S''^\ S^'^\ and S^'^\ A suitable parameterization of the 
decomposition of the Fourier components of T^^'' is 



(2.43) 



^00 


— /pj 








f Wj, 


rp{<j) 

ij 


~ ^ij fp 





j fn + -jiWtkj + Wjki) + Tij, 

where /p, /„, /p, and /^r are random function of k; w and v are transverse vectors, w • k = v • k = 0, and 
is a symmetric, traceless, transverse tensor, = fc-' = 0. The variables (/.), (v,w) and (r^j) represent the 
scalar, vector, and tensor degrees of freedom of T^^^ respectively. They are the source of the perturbation 
equations. 

The goal of the next three subsections is to express the correlators of the source components 
and iS^"'"-' in terms of these variables. These expressions, inserted in the perturbation equation ( 2.36| ), then 
allow us to compute the CMB anisotropy and dark matter power spectra numerically. 



D. Axion seeds — Scalar component 



We first consider the scalar contribution given by the four variables /, of Eq. |2.43| . Only two of these 
functions are independent, the other two are related by energy and momentum conservation. We shall use 
two linear combinations of the three scalar seed- functions fp, fy, and /tt, 



/p(k,r;)=aV-)=ro(o^)(k,,7) 



3 

2fc4 



(2.44) 
(2.45) 
(2.46) 



In the presence of seeds and in the linear perturbation approximation, the scalar component of the total 
geometric perturbations determined by the Bardeen potentials $ and 5* can be separated into a part induced 
by the seeds, and <I>s, given by 



kHs = 4ttG [fp + 3(d/a)/„)] , $s + *,s = -SnGf^, 



(2.47) 



and a part induced by the perturbations of the cosmic fluid, "^fm and $,„ . The total geometric perturbations 
are given by the sums. 



(2.48) 



The Bardeen potentials are gauge invariant and fully describe scalar perturbations of the Friedman geometry 
(for details look in HjH). 

Scalar perturbations are seeded by and ^'s. These are the standard independent variables to use as 
scalar sources in the perturbation equations. In order to simplify somewhat the computation, we use <I>s and 
fjr as our scalar seed degrees of freedom and the scalar source vector becomes 
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5(^)(k,,7) = [$,(k,7?),47rG^(k,r;)]. 



(2.49) 



The energy- momentum tensor of the axion is given by Eq. ( ^.35| ), which leads to the foUowing expressions 
for the seed- functions in terms of the axion field a: 



o'(P,f?)o-(|k- p|,?/) - p • (k- p)CT(p,?/)CT(|k - p|,7y) 



1 f d*p 

f^(k,ri) = -— J ^^-^k- (k-p)cr(p,77)cr(|k-p|,77) 
d^p 



(27r)3 



(k • p)[k • (k - p)] - -Pp ■ (k - p)J a{p, r^)a{\k - p|, ,7). 



(2.50) 
(2.51) 
(2.52) 



The first two seed- functions, fp and /„, together with Eq. (2.47), yield 



d^p 
(2^ 



-CT(p,?7)(T(|k- p|,77) - -p- (k-p)cr(p,r7)CT(|k-p|,77) 



AttG 



„a k • (k — p) . , , ,,, , , 
3 4:^a(p,r;)a(|k-pU) 



(2.53) 



The only information about the source random variables which we really need are the unequal time 
correlators between the Fourier components of the independent variables $s and /^r. These correlators can 
be written in terms of four real (since the correlators ((T(k, 77)17* (k', 77')) are real) scalar source correlation 
functions, Fn, F22, F12, and F21, which completely characterize the scalar component of the source, 

(<i>,(k, 77)$: (k', 77')) = 5(k - k')i^ii(fc, 77, r;'), 

4^G(<i>,(k, 77)/; (k', 77')) = (5(k - k')Fi2(fc, 77, 77'), 

4^G(A(k, 77)$: (k', 77')) - (5(k - k')^^2i (fc, V. ri'), 

{inG fiUiK V)f: (k', v')) = - k')F22{k, 77, rj'). 

Note that Fii{k,ri,ri) and F22{k,r],ri) are positive by definition and, since the functi ons F, are real, 
Fij{k,ri,r]') = Fji{k,ri' ,ri). In order to compute these functions we make use of Eqs. (|l5|) and ( |I53| ) 
and we exploit the stochastic average conditions of the Gaussian variables a and a (Wick's theorem). We 
first introduce three real auxiliary variables Si, S2, and E3, which depend on the power spectrum of the 
axion field, (|c(k)p), and on the solution ip of the evolution equation, Eq. (2.32), 

(a(k, 7?)a(k', 77')) = (27r)3<5(k - k')Si(fc, 77, 77'), 

(a(k, 77)a(k', rj')) = (27r)-^(5(k - k')S2(fc, 77, 77'), 

(a(k, 77)a(k', r/)) - (27r)35(k - k')i:,ik, 77, 77'), 

(a(k, 77)a(k', r/)) = (27r)3^(k - k')^3{k, t?', 77). 



The variables are given by 



^i{k,V,v') = 
^2{k,V,v') = 
^3{k,V,v') = 



ka[ri)a{ri') 

/L^WH [^(fc, ^) HivMk, 77)] [^{k, 77') - niv'Mk, rj% 

Ka[rj)a[rj' ) 



|c(k) 



ka{r])a{r]') 



'-[^{k,fj)~H{v)Mk,fj'), 



(2.54) 

(2.55) 
(2.56) 
(2.57) 



where Ti. = d/a. Notice that Si (77, 77) an d ^2(7 7, 77) are positive by definition. 

Inserting these results in Eqs. ( ^.52 ) and (2.52), and making use of Wick's theorem for the "random 
variable" c(k), we can work out a somewhat lengthy but straight forward expression for the scalar source 
functions, Fn, F22, F12, and f2i, in terms of the variables Ei, E2, and S3: 
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fc4 



/ (^{^S2(p,r;,,y')S2(|k-p|,77,r?') 



--P • (k - p) [SsCp, 77, 7/)I]3(|k - p|, r;, 77') + E3(p, 77', 77)E3(|k - p|, 77', 77) 



k.(k-p) 



fc2 



H(77)S2(P, ry, 77')S3(|k - p|, 77, r;') + 7^(r7')S2(p, r?, 77')S3(|k - p|, 77', 77) 



i(p . (k - p))2S]i(p, r;, 77')Si(|k - p|, 77, 77') + g lP-l^ P^)if P-k) ^ 
H(77)S3(p, 77', 77)Si(|k - pU, 77') + 7^(77')S3(P, V. V)Si(|k - pU', 7;) 



fc4 



k.(k-p))2S2(p,7y,77')Si(|k-p|,7;,77') 
(k • (k - p))(k • p)I]3b, ?7', ^7)S3(|k - p|, 77, 77')] } , 



i^i2(fc,77,7;') = 



(27r)3 



(k • p)(k • (k - p)) - ifc^p . (k - p)] 'Si (p, 77, 77')Si(|k - pU, 77') 
j-^ [(k • p)(k . (k - p)) - ifc2p . _ p) 



E3(p, 77', 77)S3(|k - p|, 7;', 7;) - p • (k - p)Si(p, 77, 77')Si(|k - p|, 77, 7/) 
k-(k-p). 



-I]3(p,7y',77)Si(|k-p|,7;,77') 



F2i{k,Tl,'n') = F^2{Kti' ,11). 



The scalar source correlators of the perturbation equation (2.36) can be written as a two by two positive 
and hermitian matrix, 



(5f^(k,77)5f)*(k,77')) = 



Fll(fc,77 , 77') Fl2(fc,7;,77') 
F2i(fc,7;',77) F22(fc,77 , 77') 



(2.58) 



E. Axion seeds — Vector component 

The vector contribution to the perturbation equations is seeded by the vector seed-functions Vi, Eq. ( 2.43| ), 

7;.(k,77)=i-fro(;)(k,77), (2.59) 
where is the projector operator onto the space orthogonal to k defined by 

Pij ~ k^kj, fcj — ki/k. (2.60) 

Again, the second vector seed function, w, is given by v via momentum conservations. Defining the projection 
of the vector p onto the space orthogonal to k by p^ — Pp, we obtain an expression for the vector seed- 
functions in terms of the axion field, 

j—j^pfa{p,ri)a{k-p,T]). (2.61) 

We again need the unequal time correlators between the Fourier components of the vector seed-functions 
Vi. These correlators can be written in terms of a vector source correlation function G, which completely 
characterize the vector component of the source , 

(47rG)2(7;,(k,7;)7;,(k,77')) = (<5y- - fc,%)G(fc, 77, 77'). (2.62) 
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Using Eq. ( |2.6l| ) and Eq. ( ^.54 ) this function takes the form 
G{k,i],T]') = 



2fc2 



The vector source correlators of the perturbation equation ( ^.36 ) then are 



(2.63) 



(2.64) 



F. Axion seeds — Tensor component 



The tensor contribution to the perturbation equations is seeded by the tensor seed-functions , Eq. (2.43), 

(2.65) 



r,,(k,77) = ( P^P- - tp^.P'n ) Ti:)(k,77). 



1 

2' 

This leads to an expression for the tensor seed-function in terms of the axion field. 



(k, 77) 



(27r)3 



pfpj - (l/2)(5,, - hkj){p^ f a(p,77)a(k- p,77). 



(2.66) 



which can be used to compute the unequal time correlators. These correlators can be written in terms of 
a tensor source correlation function, H, which completely characterizes the tensor component of the source 



(47rG')^(Ty (k, 77)T;,„(k, 7;')) = [diiSj,^ + SimSji - dij6i„i + k ^{Sijkikm + Simkikj - Sakjkm 

^imkikj djikikjji Sjmkiki^ ~\- k kikjkikjji\Il (^k, rf, Tj ) 

= (PaPjn. + PJlP^m " P^J Plm) H{k, 77, 77'). 

Using Eq. ( 2.66| ) and Eq. ( ^.54| ) this function takes the form 

(47rG')2 f d^p 



7J(fc, 77,77') 



4fc4 J (27r) 



The tensor source correlators of the perturbation equation, Eq. (2.36), hence are 



(k,77)>S[^)(k,7;')) = {PaPJ„,+PJlP„n~P^JPlrn)H{k,7J,^J'). 



(2.67) 



(2.68) 



(2.69) 



III. CMB ANISOTROPIES INDUCED BY AXION SEEDS 



In this section we present the CMB power spectrum obtained in our scenario. We first describe the CMB 
angular power spectrum obtained in the coherent approximation and in Subsection B we then show in detail 
that the coherent approximation is very good for axionic seeds, leading to errors of 5% or less. 
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A. CMB power spectrum — Coherent approximation 



A source is called coherent p6| , p7| if the unequal time correlation functions can be factorized or replaced 
by the product of deterministic sources, as in Eq. ( 2.42 ), 



(3.1) 

As pointed out in Subsection |II Cj , this approximation is exact only if the source evolution is linear. Then 
the different k modes do not mix and the value of the source term at a fixed k at a later time is given by 
its value at initial time multiplied by some transfer function, as in Eq. ( |2.40| ). In this situation Eq. ( |3A| ) 
becomes an equality and the model is perfectly coherent. This is not the case for our model since we know 
that, although the axion field evolves according to a linear equation, its energy-momentum tensor, which 
enters into the perturbation equations as source, does not; it is quadratic in the field a. Thus, nonlinearity 
leads to mixing of scales and to deviation from a Gaussian distribution. 

Nevertheless our situation is very similar to the large N limit of global 0{N) models in which the only 
nonlinearities also are the quadratic expressions of the energy-momentum tensor. In this case the effects 
of decoherence are very small and one finds that the full incoherent result is not very different from the 
perfectly coherent approximation 

This result motivated us to compute the CMB anisotropy in the perfectly coherent approximation. Here 
we repeat and expand on results already presented in p5|| while in the next subsection we justify them by 
discussing the full incoherent case. 




FIG. 3. Time evolution of the source function Fii{k,rj,'ri)k^ , with tilt rio- = 1.1, for different modes, k — O.OOlfee 



0. 1/Ceq , k 



Wkeq. For super-horizon modes, the correlator Fn decays like rj 



' /k" 



As soon as a mode enters the horizon the corresponding correlator decays faster due to the oscillating behavior of the 
axion field. Before crossing the horizon, the other scalar equal time correlators show the same power law behavior 
while the vector correlator G{k,ri,ri) oc rj^^^"''' /k^ and the tensor correlator H{k,ri,r]) oc ri^^'^"" (independent of k). 

In order to compute the CMB anisotropy power spectrum in the coherent approximation, we replace the 
unequal time correlation functions in Eq. ( ^.39 ) by the products 



(5(^)(k,r;)5(^)(k,r,')) - G^""^^, ^7, ^ [G^'''Hk,V,v)G^''^Hk,r,' ,^')]'^' , 
(5(^)(k,7/)5(^)(k,ry')> = H^'^'Hk^V^v') - [H^'''\k,Tj,r])H'-'''^k,r]' ,r]')]^^\ 



(3.2) 



where we have indicated the dependence of the correlators on the spectral index by a super-script. In 
Fig. we show the time behavior of one of the equal time correlators. On super-horizon scales, krj ^ 1, 
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they all display the same typical behavior, oc fc""?;-'^"^"" , which depends on the spectral index no- and on 
K, a positive power determined by dimensional arguments. On sub-horizon scales the correlators decay fast 
due to incoherent oscillations of the convolved axion field. 



We have solved Eq. (2.3(:) for the scalar, vector, and tensor components. The CMB anisotropy power 



spectrum is given by the sum of the three contributions and depends on the spectral index ria, 
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FIG. 4. The CMB anisotropy power spectrum for fluctuations induced by axion seeds with a tilt Ua = 1.1 and 
A = 0. This result is computed within the coherent approximation. We show the scalar (dot-dashed), vector (dashed) 
and tensor (dotted) contributions separately as well as their sum (solid). 

In Fig. ^ we show the scalar, vector, and tensor contributions to the resulting CMB anisotropics for an 
axion spectrum with tilt — 1.1. The "hump" at £ ~ 60 in the scalar component is due to the isocurvature 
nature of the perturbations. This is also one of the reasons why the acoustic peaks are very low, the other 
being that the vector (and tensor) component is of the same order of magnitude as the scalar one. This 
enhances, in seed models, the CMB spectrum at large scales thereby lowering the acoustic peaks at small 
scales. The result obtained is remarkably similar to the large N case studied in |^^. The main difference 
here is that, like for usual inflationary models, we dispose of a spectral index which is basically free. By 
choosing slightly bluer spectra, we can enhance the power on smaller scales. 

In Fig. I we show the sum of the scalar, vector, and tensor contributions comparing the results from 
different tilts with and without a cosmological constant. The CMB power spectra obtained can have consid- 
erable acoustic peaks at £ ~ 250 to 300, which can be raised further by adding a non- vanishing cosmological 
constant. Increasing the tilt Ua- raises the acoustic peaks and moves them to slightly smaller scales. As 
found in [ p3| , the power spectrum of the scalar component is always blue. The tensor and vector component 
counterbalance the increase of the tilt, maintaining a nearly scale invariant spectrum on large scales. The 
models can be clearly discriminated from the common inflationary spectra by their isocurvature hump and 
by the position of the first peak. A discussion on the comparison of these results with recent CMB data will 
be given in Section IV. 

We have also computed the CMB polarization for our model. The result for two different spectral indices is 
shown in Fig. |^ where we compare it with the polarization from usual inflationary models. It is interesting to 
note that our models show a characteristic "polarization hump" which is significantly smaller in inflationary 
models. The polarization "hump" is completely suppressed for topological defects due to causality and 
represents a very characteristic signature of "acausal seed models" like the one under consideration. 
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FIG. 5. The CMB anisotropy power spectrum for fluctuations induced by axion seeds. We show the sum of the 
scalar, vector, and tensor contributions for 5 difi'erent tilts, with Q.h = (solid) and Q,k ~ 0.7 (long dashed). The tilt 
is raising from bottom to top, ria = 1.1, 1.2, 1.3, 1.4, 1.5. 



B. Decoherence 



In order to estimate the accuracy of the results found in the previous subsection, we discuss here the 
decoherence of the axion seeds showing that the difference between the coherent approximation and the 
full incoherent calculation is very small. The decoherence is tested only for the scalar component of the 
spectrum, where it may lead to "smearing out" of the acoustic oscillations. Its effects on vector and tensor 
perturbations are expected to be small. 

We first introduce the property of "scaling" for the axion seeds. When working with seeds, to solve the 
problem of the enormous dynamical range^ needed to compute the Ci's from ^ = 2 to ^ = 1500, one often 
makes use of scaling properties. We call seeds scaling if their correlation function, (5(k, rj)S{'k, rj')), is scale 
free, i.e., the only dimensional parameters in Fij, G, and H are the variables 77, 77', and k themselves. As 
we have already mentioned, axion seeds are not scaling since the correlation function contains factors of the 
form [k/ki)'^ . But such a simple pre- factor can we written as 

(fc/fci)" = {kriT/ikirir 

and does not enter the costly numerical integration. Numerical calculations are reduced greatly if one can 
write the correlation function in the form 

G{Kfj,7^')=gi^,k,)W{y,r), (3.4) 
H{KV,V') = HVvv',ki)T{y,r), 

where y = ky/rjrf and r = y'r/Jr], and /, g, and h are given explicitly. The matrix Cy and the functions W 
and T are dimensionless by construction. In the following we shall call this behavior "modified scaling" . 



■^To compute the CMB and dark matter power spectra, we need to know the seed functions over a dynamical range 
of fcmax/fcmin ^ 30'000 and this for all times rji„ < g,!]' < rjo with krjin <C 1. This gives finally more than 1000 
functions of two variables which have to be known accurately over a long time interval. 
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FIG. 6. The CMB polarization power spectrum in linear scale (left) and log scale (right) for fluctuations induced by 
axion seeds and shown for 2 different tilts, with JIa = 0.7, — 1.3 (lower solid line) and na — 1.5 (upper solid line), 
are compared with the standard inflation result (dashed line) for the same cosmological parameters. Polarization 
clearly distinguishes between inflation and axion seeds, especially via the isocurvature hump. 



But even after this extraction of the explicit breaking of scaling, our source does not exactly obey "modified 
scaling" due to the radiation-matter transition . As one can see immediately from the evolution equation 
of the axions in the post-big bang phase, Eq. (2.32), the extra dimensional parameter implicitly contained 
in the unequal time correlators is ry^,, which comes from the expression for the scale factor a, Eq. (2.19). 
The radiation-matter transition introduces the new scale 77* and thereby spoils the modified scaling behavior 
of the axion seeds^. However, deep in the radiation or matter era, <C 77* or 77 ^ 77* respectively, the 
reduced correlation functions do obey scaling. In order to avoid this problem and to simplify the numerical 
calculations, we therefore compute the axion field according to the equation for the pure radiation era, i.e., 
setting a{r]) = rj. We call this the radiation approximation. This approximation affects the correlators 
and the CMB anisotropy power spectrum, especially at large angular scales, but is expected not to differ 
significantly from the correct results on the scales of the acoustic peaks, and it allows us to obtain sources 
which obey modified scaling. 

In the coherent case, where we just need the equal time correlators, the numerical requirements have not 
been very involved and we have not been pushed to the radiation approximation. But, as we shall see, the 
fully decoherent calculation will not change the results considerably and therefore an enormous numerical 
effort, which would be needed to compute the unequal time correlators without any use of scaling behavior, 
is not justified for this simple test. 

In the matter dominated era, axion seeds are amplified by quantum particle creation while in the radiation 
approximation they do not experience this amplification. Nevertheless, axions are massless particles and they 
behave like a perfect radiation fluid. Thus, their energy density decreases as 1/a^, faster than the cosmic 
fluid in a matter dominate universe, where a (x rf and p cx , than in a radiation dominated universe, 
where a oc r] and p cx a~'^. This leads one to some overestimation of the sources at 77 > 77* in the radiation 
approximation. 

In Fig. 1^ we compare the time behavior of one of the equal time correlators taking into account the 
radiation matter transition (dashed) with those obtained in the radiation approximation (solid line) for 
two different values of k. Modes that enter the horizon before matter-radiation equality, k > k^q, do not 
feel quantum particle creation; therefore, there is no difference between the full result and the radiation 



This breaking of scale invariance is also found in models with topological defects. 
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FIG. 7. Time behavior of Fii{k,r),r))k^ , 
matter-radiation equahty, ip{k = 10keq,T)) 
approximation, dashed hnes without approximation. For k > k^q there is no difference on super-horizon scales, 
while for k < keq the additional amplification experienced in the matter dominated phase is lost in the radiation 



= 1.1, for a mode which enters the horizon before 
aqytf). Solid lines show the modes in the radiation 



approximation. On sub-horizon scales, the radiation approximation decays slower than the correct result, 
behavior is found for the other correlators. 



A similar 



approximation on super-horizon scales. Inside the horizon, in the matter era the mode decays faster than 
in the radiation approximation. Modes which enter the horizon after equality, k < keq, get first amplified 
by particle creation, an effect which is missed in the radiation approximation, but then decay faster than 
in the radiation approximation. As can be seen in Fig. |8[ the slower decay has consequences on the CMB 
anisotropy power spectrum: using the radiation approximation somewhat enhances the Sachs- Wolfe plateau 
and the first peak. 

We now compute the CMB anisotropics in the full decoherent case for the radiation approximation, 
making use of modified scaling. We restrict our attention to the scalar component, where decoherence can 
be important. 



As explained in Eq. (3.4) we write the scalar correlation matrix Fij (for = 1) as 

F,,(k,rj,7j') = {rjrj'f/'a,{y,r), (3.5) 

where is only function of y and r and hence dimensionless. The matrix C,j is clearly symmetric under 
r — > 1/r as can be seen in Fig. ^. For y < 1 the sources decay like 1/y and after horizon crossing they begin 
to decay faster due to oscillations. 

The source correlation matrix Cij can now be considered as kernel of a positive hermitian operator in the 
variables x = krj = y/r and x' = krj' = yr, which can be diagonalized |]3^ , 

Cij{x,x') = ^ A„Ui„(x)wj„(a;'), (3.6) 

n 

where {wj;„} is an orthonormal series of eigenvectors (ordered according to the amplitude of the corresponding 
eigenvalues) of the operator Cij for a given weight function w. The eigenvectors and the eigenvalu es d epend 



Eq. ( |239| ) leads to 
(X,(k,r;o)X,(k,77o)) = ^ A(")xf )(k,,7o)^i"^(k,7yo). 



on the weight function w which can be chosen to optimize the speed of convergence of the sums (3.6) 
Inserting Eq. (gj) in Eq. ( |2.39| ) leads to 



(3.7) 
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FIG. 8. Scalar contribution to the CMB angular power spectrum, computed in the pure radiation approximation 
(solid line) and without approximation (dashed line), with an axion spectral index ria = 1.1 and A = 0. 



where x|"'(7yo) is the solution of Eq. ( 2.36| ) with deterministic source term 
For the scalar CMB anisotropy spectrum this gives 

N 
n=l 



(3.8) 



(3.9) 



is') 

C\ is the scalar component of the CMB anisotropy induced by the deterministic source f „ and N is the 
number of eigenvalues which have to be considered to achieve good accuracy. 
In our model we actually find it easier to diagonalize the matrix 



')\/xx', 

whose diagonal is flat for x < 0.01, exactly as in the large- and texture models studied in 
we have 

N 



^ V V 



In this case 



(3.10) 



where Vjn and A„ are the eigenvectors and the eigenvalues of the matrix Cij. 

We diagonalize the matrix Cij using the logarithmic weight function w = l/x which allows us to sample 
the range of scales of interest more evenly. In Fig. |l^ we show the eigenvectors decomposition of one of 
the scalar correlators. Note that a rather high number of eigenvectors and eigenvalues is required to reach 
a good accuracy in the approximation of the diagonal of the correlation function. Summing up = 50 
eigenvectors the convergence is guaranteed; the summed up correlation function reproduces the original to 
better than 1%. 

This is different from the large- model, where about 20 eigenvectors suffice for the same accuracy. We 
assume that this difference is due to the slower decay of the source functions. As can be seen from Fig. nu, 
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FIG. 9. The correlator Cii(j/, r) is shown. In the left panel the solid, dashed, and dotted lines respectively represent 
Cii(l X 10~^,r), Cii(l X 10~*,r), and Cii(0.03, r). In the right panel the solid, dashed, and dotted lines respectively 
represent Cii{y,l), Cii(?y,0.3), and Cii(?y,0.1). The other scalar correlators C'22 and C12 behave similarly. 



the source function is decaying from its original value to about 1% over the interval 0.1 < ki] < 10, while in 
the large- iV model this decay is achieved in the interval 0.5 < krj < A. 

We now compute the scalar contribution to the CMB anisotropics using Eq. (p.9\). The result is shown 



in Fig. 11, We note that decoherence slightly reduces the amplitude of the oscillations around the first peak 
leaving however the secondary peaks and their positions almost unaffected. Although axion perturbations are 
in principle incoherent, it is difficult to observe this from the CMB power spectrum. The effects of decoherence 
are indeed very weak and the spectrum obtained in the perfect coherent approximation reproduces the 
decoherent result within less than 5%. We hence are confident to obtain a sufficient accuracy in the perfectly 
coherent approximation which we shall apply for the rest of this paper. 



IV. COMPARISON WITH CMB ANISOTROPY DATA AND MATTER PERTURBATIONS 

In this section we compare the results found in the previous section with data discussing in particular 
the consequences of the normalization of CMB anisotropics to COBE scale and presenting the cosmological 
parameters favored by our model. In Subsection D we finally compute the dark matter power spectrum and 
we compare it with data. 



A. Normalization and the kink 



Comparing our numerical result with the CMB data we normalize our curve to the fluctuation amplitude 
observed by COBE. This provides a relation between the string and the scale of the break ki,. Since we ignore 
constant factors of order unity in the overall amplitude in our calculation, the result for the amplitude 
is not very precise, but certainly correct within a factor of about 2. For the best fit value of the tilt, 
- 1 = £ ~ 0.33, our numerical result on the COBE scale (at i - 10) is i{i + l)Ci ~ O.Sgf (ry^fc;,)-^^ 
Here gi is the dimensionless string coupling constant given by wi/mpianck where uji = fci/ai = H{rii) is the 
inverse string scale. Comparing this with the COBE normalization, i{l + 1)CiTq ~ 5225/iK^, yields 

r/,fcb = (2.1 X l^y'glY'^ (4.1) 

For example, if the string scale is lO^^GeV, so that gi ^ 0.1, we get fcb ^ /i^/(2kpc), where we have inserted 
77* ^ 20/i^Mpc. An interesting constraint comes from the fact that the break in the spectrum should be on a 
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FIG. 10. The sum of the first few eigenfunctions of Cii(a::,a;') is shown for a weight function w = 1/x. The first 
(dot dashed), first and second (short dashed), first ten (long dashed), and fifty (solid) eigenfunctions are summed up. 
The open circles represent the full correlation function. Here we only show the equal time diagonal of the correlation 
matrix but the same convergence behavior is found in the Ce power spectrum which is sensitive to the full correlation 
matrix. 



scale which is smaller than the scale represented by the first acoustic peak in order not to reduce the latter. 
Since 77, corresponds to the horizon scale at equahty, this requires T^^fcf, > 1 or a;i(ai) = Hi > 0.02mpianck- 
Together with Hi < O.lmpianck, this brackets the string scale just in the bull park where it is expected for 
very different theoretical reasons. 

The length-scale/energy-scale corresponding to the break at the time 77^, during the pre-big bang phase, 
when the expansion law is supposed to change, is given by 



\h\-\VbHm)/a,^\Vb\^lO~''-\Vb\ 

aim) 



-1/4 



lO"-^-^ - 6 X lO'^^Ti - 3GeV"\ (4.2) 



where we have used t?;, ^ 77* ~ 20Mpc and 771 ~ O.lcm. The energy scale obtained in this way is uncertain 
with a factor of about 10. 

In Fig. |l^we show the dependence of the CMB anisotropy spectrum on the position of the break. Typically, 
the break lowers the second and subsequent acoustic peaks while does not substantially affect the first peak. 



B. Cosmological parameters 

In the last two years, a peak in the CMB power spectrum at ^ ^ 200 as been detected by several 
different experiments, most notably TOC098 |l[, B97 [||, B98 [|, and MAXIMA-1 |]. Among them, the 
BOOMERanG-98 power spectrum reported the best and at the same time most conservative detection, 
although coming from only 5% of their overall dataset. The position, amplitude and shape of the peak can be 
fitted by the power spectra expected in the simplest inflationary scenario based on adiabatic perturbations in 
a spatially flat universe Therefore, this peak represents the biggest challenge for the model presented 

here. 

We want to investigate whether a suitable choice of cosmological parameters can bring our model in 
agreement with the above mentioned data. This question is also very important in view of the usual 
"determination of the cosmological parameters" from CMB anisotropics, in the sense that it shows how the 
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FIG. 11. The scalar contribution to the power spectrum is shown for a primordial spectral index Wo- — 1.1. 
From bottom to top, the sohd lines show the contributions of the sum of the first ten, first twenty, first thirty and 
first forty C^^"'''s. The thick solid line represents the full eigenvectors summation (up to N — 50) to be compared to 
the perfect coherent approximation, shown by the dashed line. The decoherence does not significantly wash out the 
acoustic peak and the oscillations. 



results can change when assuming a different model of structure formation. In other words the so called 
"measurements" of cosmological parameters from CMB anisotropics are strongly model dependent! 
The peak position is determined mainly by the angular diameter distance parameter 




R- 



Here fix — 1 — — is the curvature parameter and 

{sinh y (open) 
y (flat) 
sin y (closed) 

depends on the geometry of the universe. The variable y is the following integral: 



(4.3) 



y 



(4.4) 



(4.5) 



As pointed out in | p9| , the condition R — constant identifies curves in the Qm — plane, with nearly 
degenerate Cg spectra, providing that the baryon density parameter Jlbaryon is kept constant. 

In Fig. 13 we plot likelihood contours, obtained as follows: we rescale the string cosmology power spectra 
plotted in Fig. ^ both in amplitude A (in COBE units) and position R. We compare the resulting spectra 
with the BOOMERanG and MAXIMA-1 data in the region up to ^ < 400 by a simple x^-At. We find that 
the 68% confidence limit for R marginalized over A is 1.50 < R< 1.63 with R = 1.57 as best fit (see Fig. [l3| ). 

In Fig. |lj the confidence levels on R are translated to confidence levels in the — ilm plane which are 
then combined with the current SNla results |^^. It is clear from this figure that the model can be brought 
in reasonable agreement with observations only if the universe is closed. The deviation from flatness becomes 
less and less important towards fim 0, where all the R = const lines converge at n\ = 1. While the region 
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FIG. 12. The influence of the break position on the CMB power spectrum. The top solid line is the spectrum 
without break. The dashed lines from top to bottom represent a spectrum with break at = 3/r;*,2/r;, and l/r;« 
respectively. 



with fl„i > 1 can be safely excluded from different cosmological observations, a moderately closed universe 
with Oa ~ 0.85 and fi^ ^ 0.4 is compatible with SNla results and also with estimates for from cluster 
abundance and X-ray data (see e.g. 0). 

As we have seen, the position of the first acoustic peak can be adjusted by choosing JIa and ilm so that the 
resulting universe is marginally closed. Nonetheless, the width of the peak, compressed by the increase of R, 
is still not in very good agreement with the data, as well as the isocurvature hump. The resulting normalized 
is about ~ 1.8 for the best-fit, which "excludes" the model at 70% confidence. One has however to keep 
in mind that the C'i's are not Gaussian and therefore the probability for our model to lead to the measured 
CMB anisotropics is even somewhat higher than 30%. In Fig. |l^ two theoretical CMB spectra normalized to 
the COBE data are shown together with the MAXIMA and BOOMERanG98 data. We did not optimize on 
the axion spectrum, or the baryon density parameter, but we chose n^r = 1.33, f2„i — 0.4, and r^baryon = 0.05. 

Playing with the break-scale kh we can in principle lower the second peak leaving the first one almost 
unchanged. Nevertheless, the position of the second peak is different from the one indicated by inflationary 
models and the data. Inter-peak distance is therefore a better estimator of the validity of a model. Clearly 
more and better data around the isocurvature hump region, i.e. £ ^ 100, is needed to decide definitely 
whether the model is ruled out. This will most probably be achieved with the MAP satellite |Q planned 
for launch in 2001. 



C. Polarization 



The polarization spectrum distinguishes easily between adiabatic inflation and the axion seed model (see 
Fig. |l6| ) . The preferred closed universe for axion seeds translates into a smaller distance between polarization 
peaks. As the physical distance between peaks depends only on the sound speed, which is only slightly 
dependent on f^baryon^^, a quantity which is already tightly constrained by nucleosynthesis, the A£ on which 
this distance projects is mainly determined by spatial curvature, fix (it depends also somewhat on Q\ as 
can be seen from Eq. (4.5)), and is independent on the model for the initial fluctuations. 
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FIG. 13. Confidence levels (68%, 95%, and 99%) for the rescaling factor R and the amplitude in COBE units A, 
from the recent BOOMERanG and MAXIMA-1 observations. 

D. The dark matter power spectrum 

The computation of the dark matter power spectrum had aheady been performed in p5| ] where a consid- 
erable deviation from the data was found. In this work we repeat this computation taking into account the 
preferred values of the axion spectral index and of the matter energy density found from CMB data, and 
we introduce the break in the axion spectrum discussed above. With this additional input it is possible to 
establish reasonable agreement between the data and the dark matter power spectrum (see Fig. p^ ). 

Since the computation of the theoretical matter power spectrum for a closed universe is relatively involved 
and since, for the purpose of comparing the theoretical spectrum with observations, we are interested in scales 
much below the curvature scale, we have computed it for a flat universe, with matter and a cosmological 
constant, assuming that the contribution from curvature is negligible on the scales under consideration. 
Indeed, what really plays a role for the matter power spectrum is the matter content, r2„i, which fixes the 
time of equality between matter and radiation, determines when structures can start growing, and fixes the 
position of the bend in the power spectrum. 

In Fig. we present the theoretical dark matter power spectra together with the data as compiled by 
Peacock and Dodds |Q. Depending on the scale of the break in the axion spectrum, Ab = 1/fcb, our model 
can be compatible with data for different values of Q.m in the range 0.2 < Vtm < 0.4. The role of the break is 
the following: if Ab is small we subtract power only from small scales and we are able to reproduce a power 
spectrum in good agreement with data even if is relatively high. However, if we do not introduce any 
break in the axion spectrum we find too much power on small scales and our theoretical dark matter power 
spectrum is incompatible with data (compare our present result with those found in p5[|). 

The root mean square mass fluctuation within a ball of radius Sh~^ Mpc for the model with ric = 1.33, 
fcf, = 3/77*, and Q.m = 0.25 and for the model with = 1.33, kf, = 1/f?*, and fi^ = 0.4 are erg = 0.85 and 
CTs = 0.74 respectively. Analysis of the abundance of galaxy clusters suggests as ^ 0.5r2~''-^ Q. 

E. Conclusions 

We have shown that it is possible to choose cosmological parameters which bring our model in reasonable 
agreement with the present CMB anisotropy measurements, which is however less favorable than the striking 
fit of simple flat adiabatic inflationary models. This is our main result. 

Even if our model will turn out to disagree with better data, we believe that we learn the important 
lesson that cosmological parameters obtained from CMB anisotropics are strongly model dependent, a point 
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FIG. 14. The 68%, 95%, and 99% confidence levels for the cosmological parameters f^A and Qm, from the peak 
position detected by BOOMERanG and MAXlMA-1 for the model presented in this paper (dashed). The solid 
contours are obtained including the supernovae data. 



which is swept under the carpet by the vast majority of the circulating "parameter-fitting" literature. We 
believe that it is very important in the future to concentrate on model independent quantities, like inter-peak 
distances, to determine cosmological parameters. 



V. GRAVITATIONAL WAVES 



Gravitational waves represent one of the most powerful tools to investigate the early history of the universe. 
They decouple at a temperature comparable to the string scale which makes them an important window 
for cosmological phenomena related to the string theory domain. In this section we show that axions can 
contribute substantially to the production of the gravitational wave background in the pre-big bang model, 
acting as a source in the tensor perturbation equation. This leads to a spectrum which is different from the 
standard gravitational wave background of string cosmology based on the "direct mechanism" of graviton 
production by amplification of quantum vacuum fluctuation. This new "indirect mechanism" leads to a flat 
spectrum and can easily be distinguished from the direct one. Indeed, as we shall see, the axion induced 
gravity wave background dominates the "direct background" at small frequencies and represents an important 
observational constraint for string cosmology. 



A. Direct production — Amplification of vacuum fluctuations 



So far, amplification of quantum vacuum fluctuations have been considered as the principal mechanism 
for the production of gravitational waves during the pre-big bang phase |45|-|47|]. During the dilaton era, 
before the big bang, when the scale factor evolves according to E q. ( |2.9| ) , the Fourier modes of metric tensor 
perturbations satisfy an evolution equation similar to Eq. ( 2.13| ), namely 



€ + [k' 0, (5.1) 
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FIG. 15. Two theoretical CMB anisotropy spectra normalized to the COBE data, with Qa = 0.4 and axionic 
spectral index = 1.33, are compared with the MAXlMA-1 and BOOMERanG98 data. From left to right, our 
model has a break a.t kb = S/rjt and fc^ = l/rj, respectively. Lowering ki, we subtract power on small scale and we 
can lower the second peak leaving the first one almost unchanged. 



where ar = ae "^/^ is the pump field of gravity waves and V'k canonical variable for tensor modes 

of the metric. For the isotropic case discussed in this work, one find that ar cx; jryl^^^ independently on 
the evolution and number of dimensions during the pre-big bang phase. After proper normalization to the 
incoming vacuum, this yields the solution 

V-k = (fc^), v<-m- (5.2) 



After the big bang, in the radiation dominated era, rj > r/i, the solutions of Eq. (5.1) are simple plane 



waves. From the m atching conditions between these two regimes, applying the same procedure as discussed 



in Subsection 11 B for the axion field, one obtains the following spectrum of gravitational waves. 



^,-JTT^(-X-sl(-)'^,, (5.3) 



which is a tilted spectrum, cx uj'^, normalized to at the strin g sc ale. One actually supposes that, at a string 
epoch < —771, the dilaton- vacuum regime behavior of Eq. ( ^.S| ) breaks down and the universe undergoes a 
De Sitter expansion with linearly growing dilaton, which lasts until the beginning of the radiation dominated 
era 771. This phase leads to a nearly fiat gravitational wave spectrum at very small scales. The normalization 
of the spectrum to the string coupling gi can then be performed at a lower frequency, < 1^1, leading 
to a somewhat higher density of directly produced gravitons than the one discussed here. This is very 
important in order to make the direct background observable and still compatible with nucleosynthesis. 
(See [^sf and references therein for more details.) A more detailed discussion on the important signatures 
and observational consequences of this direct production of gravitational waves can be found in |^ and 
references therein. 



B. Indirect production — Axion source 



Let us discuss now the production of a stochastic gravitational wave background generated by the presence 
of axion seeds. This indirect background will be superimposed to the direct one discussed above and will 
dominate the total spectrum at large scales. These two production mechanisms are fundamentally different. 
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FIG. 16. The CMB polarization spectrum of our model (solid line) for the best fit parameters is compared with 
the inflationary CMB polarization spectrum in a critical universe with Qa = 0.7. The fact that in our model the 
universe is closed is visible in the smaller distances between successive peaks. 



While the direct production of gravitons takes place during the pre-big bang phase and is due to the 
amplification of vacuum fluctuations, the indirect production is sourced by the axions during the post-big 
bang era. 

The creation, propagation, and damping of gravitational waves in a Friedman background are described 
by the tensor perturbation equation (see e.g. ||50[|), 

hij + 2—hij — Ahij = WnGa^Tij , (5-4) 

where tensor perturbations in the metric are parameterized by the traceless, divergence-free, symmetric 
tensor field hij, 

ffM- = 9^.. + K = Q = ^uh";, (5.5) 

whic h is a gauge invariant variable. As before a dot denotes the derivative with respect to conformal time. 
Eq. (5.4) is a wave equation with source term Ty. 

The tensor field hij is usually decomposed into two polarization states as 

%(x,r;) = /iX(x,,7)e=^.(x) + /.+ (x, ,/)e+(x), (5.6) 

where e^j = eje^ — efe| and e^- = ejej + e?ej are the polarization tensor fields and (e^,e^,e'^) is a local 
orthonormal basis (the wave is propagating in the 63 direction). 

The energy density of gravitational waves is given by the 00-component of the energy momentum tensor 
of the wave. This can be defined as a space-average over several oscillations, 

_ (hijh^^) _ {ill) + {hi) 

We decompose hx and /i+ in Fourier modes, 

^e^^-hxiKr^), A=x,+; (5.8) 

therefore 
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FIG. 17. The linear dark matter power spectra for fluctuations induced by axion seeds with spectral index Ua = 1.33 
and a break in the spectrum at (a) fci, = 3/?7, and (b) kh — 1/rj,, for a flat universe with flm = 0.4 (dotted), fim = 0.3 
(dot dashed), and ~ 0.25 (dashed) are compared with data. We assume an IRAS galaxies bias of bj = fl^'''^ . 
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'/lA(k,7?). 



The spatial average then becomes 



{hi) - / -^#^e-(''+>^')(/,,(k,^)/,,(k',,,)), 



(27r)3 (27r)3 

and we can use the stochastic average condition 

(/iA(k)/iA'(k')) = {2^f5\\i~-k!)5xx'\hx{\i)\\ 
which yields, under the hypothesis of statistical isotropy, 

1 



Pg 



dkk^\hx{k,'q)Y 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



We now compute the spectrum \hx{k,ri)\'^ in the coherent approximation. For this we introduce the 
deterministic source function n(fc, 77) defined by 



4^Ga2n(fc,77) = /hTMm?), 



(5.13) 



(for the function H, see Eq. (2.68)). The polarization tensors satisfy e^j^y — 2(5^ and we can hence rewrite 



Eq. (5^) in momentum space as 



hx + 2-hx + k^hx = 87rGa2n. 



(5.14) 



The factor 1/2 comes from the fact that 11 sources both modes x and + of hij and, assuming again 
statistical isotropy, each mode is sourced with the same strength. Since we want to compute a gravitational 
wave spectrum we only consider modes which enter the horizon in the radiation dominated era, k-q^, > 1 and 
a/ a ~ 1/77, the other modes being uninteresting (too large wavelength) for possible observations. Therefore 
we also consider modes far from COBE scale, k ^ kb, and we can comfortably assume a flat axion spectral 
index, n„ = \. We then write Eq. (5.14) as 



h'i + -h'^+ hx 



'A ^ -"-A 

X 



J{k)/^ X <1 (active source) 
x>l (dead source). 



(5.15) 
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where the conformal time derivative has been replaced by the derivative with respect to a; = kr]. In this 
equation we assume that the axion source can be approximated by a power law behavior outside the horizon 
which is of the form 



SnGa^nik, x) = x-^/^k^ f{k), f{k) ~ STiglk-^/"^, 



(5.16) 



and can be considered negligible inside the horizon where the correlators decay quickly. 

The homogeneous solutions to this equation are the spherical Bessel functions of index zero, ^0(2;) and 
2/0(2;). In the regime, a; < 1, the solutions can be found with the Wronskian method, which yields 



hx{k,x) = f{k)[ci{x)jQ{x) + C2{x)yo{x)], x < 1, 



where 



ci( 



(x) = / dxx^^^cosx, C2{x) — / dxx^^"^ 



' smx, 



while in the second regime, a; > 1, they are a linear combination of the homogeneous solutions, 

hx{k, x) = A{k)jQ{x) + B{k)yo{x) x>l. 
By matching Eqs. ( |5.17D and ( |5.19| ) at a; = 1 we find 

hx{k, 7,) = /(fc)[ci(l)jo(M + c2{l)yo{kv)], A{k) = /(fc)ci(l), B{k) = /(fc)c2(l) 
which yields, for x ^ 1, h\ ^ kh\ ^ f{k)/r], and thus 



\hx{k,v)\ 



2^ (87r)V^-3^ 



Using Eq. ( 5.12| ) we hence find 



Pa 



dk 



dpg 



dlogk irGa^rf 



(5.17) 

(5.18) 

(5.19) 
(5.20) 

(5.21) 
(5.22) 



which corresponds to a flat spectrum of gravitational waves. 

On the other hand, at early time the radiation energy density, dominates the Friedman equation which 
becomes 



SttG 



(5.23) 



With d/a ~ l/?7 we can write the gravitational wave background spectrum produced by the axion field as 



Pi 



(5.24) 



C. Observational consequences 



In the previous subsection we derived the spectrum of gravitational waves induced by axion seeds and 
we found that it is flat on scales much smaller than the COBE scale and normalized such as to lead to the 
correct amplitude of fluctuations in the CMB anisotropies. 

Its normalization depends on the fundamental ratio between the string and Planck mass which is usually 
taken to be of the order of gi ^ 0.1 ~ 0.01 |^^. The energy density of induced gravity waves is proportional 
to the forth power of gi like the CMB anisotropy spectrum. Since the COBE normalization also depends 
on fcb (see Eq. (4.1)), which plays no role for the gravity wave spectrum on the scales considered here, gi 
alone is still allowed to vary in the range cited above even though Eq. (4.1) provides a precise constraint 
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FIG. 18. The gravitational wave spectrum in pre-big bang model for a value of gi — 0.03. The directly produced 
background (dashed line) has been normalized to string mass at a frequency of = 500Hz. This frequency corre- 
sponds to the time ris — —l/{2nfsa{rjo)) for the transition between the dilaton-dominated regime and the De Sitter 
phase in the pre-big bang era (see above). The solid line represents the sum of the direct and indirect production of 
gravitational wave background. The analysis has been limited to frequencies / S> ft, where denotes the frequency 
corresponding to the scale of the break. 

for a combination of gi and kb- Using the previous values for gi we find a flat spectrum of gravity waves 
with h^rig ~ 4 x_(10~^ 10~^^), a range which, most probably, will be reached by the third generation 
interferometers ||51| . This renders the indirect gravity wave background an important observable of string 
cosmology. Note also that in the case of its detection it would provide a direct measurement of the string 
scale! 

At present the most relevant observational bound for a gravity wave background comes from pulsars. 
In particular, the timing of the milli-second binary pulsar implies a limit on any stochastic gravity wave 
background of h^ng{atf — 4.4 x lO^^Hz) < 1 x 10^^ (at 95 % c.l.) ||5^, which transforms in our case into a 
limit on gi < 0.07 in this model. 

The direct gravitational wave background has a blue spectrum and therefore dominates the indirect back- 
ground on small scales, as shown in Fig. 1^. The crossover frequency ujc between the two regimes is 
determined by gi and the normalization frequency lUs discussed above. 



2 /3 

= gi ujs 



(5.25) 



This crossover may actually, depending on the unknown value lUs, fall into the range of frequencies at which 
interferometers will be operating. 

Finally, we would like to point out that, like the CMB anisotropics of this model, the indirect gravity wave 
background considered here is not Gaussian, which can lead to interesting observational consequences. 



Sensitivity curves for LISA and LIGO are based on |5iy52[ and references therein. We acknowledge Carlo Ungarelli. 
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VI. CONCLUSIONS 



We have investigated the consequences of axion seeds which naturally occur in the context of string 
cosmology. We found that these seeds may induce the observed large scale structure and CMB anisotropics 
in the universe provided that there is a break in the primordial axion power spectrum which from slightly 
blue on very large scales turns to a flat spectrum on scales smaller than the break, k > kb. Such a break 
appears if the expansion law undergoes a transition during the pre-big bang phase. For the scenario to agree 
with observations the break must occur at rjb — 0.3?7*, which corresponds to an energy scale of the order 
of several GeV. 

The axion seed model leads to isocurvature fluctuations with important contributions from vectors (about 
50%) and tensors (about 15%) on large scales. The first acoustic peak in the CMB anisotropy power spectrum 
is around £ ~ 300 for a flat model, il = fl\ + ilm = 1. To reproduce observations the universe has to be closed 
with parameters, JIa ~ 0.85 and ilm ~ 0.4. This parameter choice is also in agreement with supernovae and 
cluster data. Even though our model leads to a larger when fit to the CMB data it cannot be excluded 
by the presently available data. However, the "isocurvature hump" at ^ ~ 40 and the reduction not only 
of the second but also of the third acoustic peaks are signatures which clearly distinguish the model from 
standard inflationary scenarios. Furthermore the CMB polarization spectrum significantly differs from the 
inflationary result. 

We have also studied gravitational waves which are generated during the post-big bang phase by the tensor 
type anisotropic stresses in the energy-momentum tensor of the axion field. We found that they lead to a 
flat observable background of gravity waves which can give stringent constraints on the model if detected by 
the planned LIGO-III and LISA observatories. 

As the model studied is very predictive let us finally mention that its failure to reproduce observational 
data, which is hinted by present CMB anisotropy measurements and might be reinforced by future more 
accurate data, does not by itself rule out string cosmology. An additional important hypothesis of the model 
is that non-gravitational interactions of the axion field with the dark matter may be neglected and the axion 
plays the role of a "seed" . If this hypothesis is relaxed, the axions may interact with radiation and dark 
matter and even lead to a standard adiabatic fluctuation spectrum. This idea deserves further study, but 
most probably the non-Gaussian character of the perturbations also survives in such a scenario. 
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